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INTEGRABLE SYSTEMS, OBTAINED BY POINT 
FUSION FROM RATIONAL AND ELLIPTIC GAUDIN 

SYSTEMS 

Chernyakov Yu. ^ 



, Using the procedure of the marked point fusion, there are obtained integrable 

' systems with poles in the matrix of the Lax operator order higher than one, considered 

i Hamihonians, symplectic structure and symmetries of these systems. Also, taking 

^ I the Inozemtsev Limit procedure it was found the Toda-like system having nontrivial 

^ ' commutative relations between the phase space variables. 
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1 Introduction 



> 

■ Rational and Elliptic Gaudin systems In the present work there were 
, used Hamiltonians, Lax operators and symplectic structure of the rational and 

elliptic Gaudin systems. These systems were studied in the works [1] — [4]. 
Let us give their short description. 
. Let us begin from rational Gaudin system and consider n marked points 

Xa,a = l,n on CP^ and assign a coadjoint orbit of the group SL{N,C) to 
each point. The coordinates p^^ are functions on the corresponding orbits. 
The Lie-Poisson brackets of (spins) have the following form (Appendix 
QJ: B.): 

{pi^,Pb} = Sab{6''p'J-S'^pii) (1) 

On the phase space, which is the direct product of the orbits factorized by 
^ ' SL(N, C) group with coordinate independent elements on CP^ , it is possible 

■ to define the integrable system. The Hamiltonians of this system commuting 
with respect to the Lie-Poisson brackets (1) have the following form: 

rr < PaPb > 

Ha = l^ , (2) 

where <> means the trace . Other Hamiltonians are also expressible in terms 
of the traces of the spin products and differences of coordinates at the marked 
points in certain degrees. 

Let us describe now the elliptic case: instead of CP^ let us consider the 
elliptic curve. The integrable system is determined on the symplectic factor- 
space R ([4]) which has the following form: 

R = {u,v)x(Y[Oa//H), (3) 

a 

where Oa - the coadjoint orbit of the SL{N, C) group representation, 
H - Cart an subgroup, 
u,v & H - Cartan subalgebra. 
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The commutation relations between v,u and p'^^ in addition to those be- 
tween p^J have the following form: 



{v\u^} = S'^, K,pjf} = 0, {n^pjf■} = (4) 

Well-known Hamiltonian of the elliptic Gaudin system describes a physical 
system consisting of n interacting particles: 

H = Y^l{vr + J2<p'>-E2{u% (5) 
where Vi - momenta, 

yij = ^* _ _ difference of particle coordinates, 

E2{u^^) = E2{u'^,t) - interaction potential of particles, represented by the 
elliptic Eisenstien functions. This function is connected with p-Weirstrass 
function: E2{z,t) = p{z,t) + 2(^(^), where z - coordinate on the torus. Func- 
tion £"2(2;, r) is expressed in terms of ^-function ( Appendix L. ): E2{z,t) = 

Note that in the previously mentioned works there are studied systems 
which are the examples of the Hitchin systems - intcgrable systems on the 
cotangent bundles to the moduli spaces of holomorphic bundles over Riemann 
curves. These systems may be obtained by the symplectic reduction ([4], [5]) 
from the input phase spaces. Using the Lax operator we get the important 
way to describe the integrable systems. So, in the case of rational Gaudin 
system the Lax operator takes the following form: 

^<^' = ?(F^- (6) 

where z - coordinate on CP^. 

The matrix of the Lax operator in the elliptic case takes the following 
form: 

diagonal part - 

V' = v' + Y.p:-E,{z-Xa) (7) 



non-diagonal part - 

a 

where: 



{z - Xa) - (Z- Xa) 



T — T 



ip{u'^,Z-Xa), (8) 



Q{u'-3)Q{z - Xa) 

In the case of the Hitchin system the matrix of the Lax operator coincides 
with the solution of the moment map equation appearing in the symplectic 
reduction procedure. The Hamiltonians of the intcgrable system come from 
the appropriate basis expansion. The Lax operators studied in the above- 
mentioned works admit the first order poles. 

In the present work we try to answer the naturally arising questions: how 
do the systems with the higher poles in the Lax operator look like, what kind 
of symplectic form do these systems have and what are the moduli spaces of 
these systems. 



Results and methods being used to obtain the systems We use 

the procedure of the point fusion to obtain the integrable systems from rational 
and ehiptic Gaudin systems. The main idea of this method consists of finding 
such decomposition of the variables p, which gives us the pole order in the 
matrix of Lax operator being higher than one at some marked point, for 
example Xa, after taking the limit — = e — > 0, where e is the parameter 
of the point fusion( we put here the coefficient at e equal to 1). It implies the 
existence of Hamiltonian and Casimir numbers in this new system being the 
same as in the initial one. This situation is realized by passing to the new 
variables using grading decomposition: 

Pa = aoP° + aip^e'^ + ... + UrP^e'"^, (10) 

where a,. - some variables. 

The Hamiltonians in the case of two point (chosen from n marked point) 
fusion [xb — > Xa) have the following form: 



fiP =< ipy > +2 Y. = t 



< p^Pc > < P^Pc > 



(X(j Xq) (Xa Xc) 



jj2,l ^ 2 < P^Pc > _|_ < P^Pc > _^ < PcPd > 

{Xc — Xa) {Xc — Xa) ^ ^ (Xc — X^ 

(11) 

where <> means the trace. The commutation relations in this case are 
determined in chapter 2, formula (22). The Lax operators of these systems 
were described in work [6]. Note, that in the case of r point (chosen from n 
marked point) fusion the Lax operator has the following form: 

^ ^ p'' ■ kr _^ _^ ■ h _^ P° ■ fcp ^ Pc 

(Z-Xa)'- (Z-Xa)"^ {z - Xa) ^^^^{z - XcY ^ ' 

where - some coefficients depending on the parameters of decomposition 
p*", the parameters of the point fusion and some conditions arising from the 
requirement of the absence of singularities in the final expression. There is a 
correspondence between p^ and upper triangular matrices which are isomor- 
phic to polynomial algebra of e^^ variable with coefficients p^ . Commutation 
relations between p^ were found in chapter 2 (table (28)). Note, that special 
case of the Lax operator (12) was obtained in the work [7]. 

In the elliptic case the Hamiltonians of the two point fusion have the 
following form(p G sZ(2,C)): 

^2,2 ^ _^^0)i2(^0)2i + E,{u) ■ [-{py^p^f^ + {py\p')^''\+ 

+{p')'\p'?^ ■ [e[{u) + El{u)] - 2v{py' 
H^'' = v' + {{p')'r ■ il^iu) + {E[{u))'] + {p^W • E[{u)+ ^^3^ 

+[{pyW - iPrW] ■ [Ei{u)E[in) + lE'lin)]- 
-{P')''{P')'' ■ [{E[{u))' + Eliu)E[{u)+E,{u)E';{u) + ^^'(u)]. 



where u,v EH belong to Cartan subalgebra, 

Ei{u^^) = Ei{v}^ ,t) are the elliptic Eisenstien functions ( Appendix A. ) de- 
termined on the complex torus with modulus r. This function is connected 



with C-Weirstrass function: Ei{z,t) = C{z,t) +2C{^)z. Ei{z,t) is expressed 
in terms of ^-function: Ei{z,t) = c?^ ln(^(z, r)). 

In the elhptic case the matrix of the Lax operator of r point fusion has 
the fohowing form: 

diagonal part - 

V' = v' + {pY ■ h ■ E[{z - Xa) + ... + {p'f ■ ~kr ■ e[''\z - x„) (14) 
nondiagonal part - 

= (jPyi . k, . ip{u'^,Z -Xa) + ... + (pT ■ K ■ ^^'^ {u'^, Z - Xa), (15) 

where kr is similar to kr, commutation relations are the same as in the 
rational case. The Hamiltonians expressed by the variables invariant with 
respect to the gauge fixing are found in chapter 3, formula (50), and commu- 
tation relations - formula (53), (54). 

The symplectic structure in the two point fusion case (rational and elliptic 
one) is described in chapter 2 

The procedure called Inozemtsev Limit ([8]) is the method permits to ob- 
tain Toda-like system ( the system with exponential type of the interaction 
) from the elliptic Gaudin system ( the system with elliptic type of the in- 
teraction ) by the combination of the trigonometric limit, the infinite shift of 
the particle coordinates and the renormalization of the interaction constant. 
The technical details are given in [8] , [9] . The Hamiltonians of the Toda-like 
system , obtained by Inozemtsev limit from the elliptic Gaudin system in the 
case of the two point fusion and commutation relations were found in chapter 
4, formulas (86), (84). Note, that the other Toda-like systems having nontriv- 
ial commutation relations between the phase space variables were studied in 
the work [10]. 



2 Systems arising from the rational Gaudin 

system 

Let us begin from the rational Gaudin system. We will consider the two 
point fusion of n marked points (n > 4) and fulfil the following coordinate 
transformation and decomposition of the variables p: 

X}) — Xa ~\~ ^ 

(16) 

Pa = oiQp^ + aip^e S pb = !3op^ + I3ip^£ \ 

where ar, (3r are some parameters. Putting (16) in the matrix of Lax operator 
(6) and taking the limit e 0, wc get: 

(Z-Xa) (z-Xa-e) ~ ^c) ' ^ ' 

During the calculation we need to put some additional conditions on param- 
eters. The reason for their appearance is the requirement of the absence of 
singularities coming from degrees of if £ — > 0. In the present case we have 
one condition ai + Pi = and, setting ckq = 1, ai = — 1, /?o = 0, /?i = 1, 
we get: 

T _ P^ \^ Pc 

-'-'new — 7 Tn + 7 7 + / , 7 7 (18) 

{z-XaY {z-Xa) ^^f, (^ " ^^c) ^ ^ 



The decomposition < L^^^ > with respect to {z — Xg) has the following 
form : 

< L^^ >= hI'^ ■ {z - Xa)-^ + H^'^ ■ {z - Xa)-^ + H^'^ ■ {z - Xa)-^ + H^/ ■ (; 

+ J2 {H','' . {z - + H'/ ■ {z - x^r') 

c^a,b 



(19) 



where: 



H 



2,4 



2,1 



E 



< p^Pc > < P^Pc > 



Xc) {Xa Xc)'^ 



H. 



2,2 



(20) 



jj2,l ^ 2 ^ ^^^^ ^ + < -P^^g > 
(^^c SJq)^ (Xc a^a) 



+ E 



< PcPd > 
{Xc - Xd) ' 



where H^''^ , H^'^ , H^'^ are Casimir operators. 

The condition J2pa = 0, which is the moment constraint corresponding 

a 

to the action of the remnant gauge group, has to change because after the 
transformation of Pa and Ph, we get: 



0, 



(21) 



The numbers of Hamiltonians and Casimir operators do not change: instead 
of two Hamiltonians and Casimir operators of the rational Gaudin system 
related to the points Xa,Xb, there appeared and being 

Hamiltonians and Casimir operators correspondingly. It is possible to calcu- 
late the following commutation relations between the new variables: 



{{P^T, {P^f^} = S^\p^f^ - 5^^{p^f\ 
{{P^T, {P^f^} = 5^{p^f^ - 5^^{p^f\ 
It may be schematically presented as: 



(22) 



(p\/)-o, (p^:p^)^p^ {p\p 







^0 ^\ 



P 



(23) 



Using the elements p^,p^, it is possible to compose matrixes of the form (in 
the case r = 2): 



P 




P 



(24) 



, which generate the parabolic algebra. Taking into account the position of p 
in the matrix P (in the case when p'^ G sl{2, C)), it is possible to show that 
the commutation relations between p^ have the following form (Appendix B.): 

{{py\ {p^f} = 6'\p'^)'''^+^ - (5^^(pi)*''+^, 
{(p°)'^ ip")''^} = ^[<^''((p°)'='' + (p°)'=+^'^+^) - s'^^upY + {py+^^'+'')], 

(25) 



where N = 2. The loop decomposition of the matrix P in the case of two 
point fusion gives polynomials of the form oi +p^e~^ . The isomorphism can 
be established by comparing commutation relations for the matrix elements 
and for polynomials. Considering the transformation (16) for x and p in the 
case of fusion of the first r (for definiteness) of n marked points, we get the 
generalized expression for Lnew , not depending on the sequence of those points 
of fusion but depending (up to the coefficients) on the number of points: 

L„„,n = r- + ... + 7 + 7 r + > , r (26) 



X, 



)2 {Z-Xa) 



Xc) 



The matrixes P, which generate the parabolic algebra, have the following 
form: 



P = 



p' 


p" 





p' 









„n-2 



P 



P 



P 



P 



P 



P 
n-2 



„n-l 



P 
n-1 







p 



p 



p 

n-2 



(27) 



Schematically, it is possible to present the commutation relations in the form 
of the table: 



p 


pO 


pi 


p' 




pu-l 


pn 




p» 


pi 


p' 




pTl-l 


p"" 


p^ 


P' 


') 
P~ 






P" 





p2 


p2 


pS 


4 

p 


























pU-L 


pU 













pn 


pU 

















(28) 



The intersection of the column and the row is the result of the commutation 



Symplectic form in the case of the point fusion Following works 
[2] — [4] let us consider the symplectic form of the Hitchin system in the case 
of n marked points on CP^, selecting two points for the next fusion: 

n= f {5A/\6A) + J2{Pa,9a^Sga) = 

a 

= / [{5A ^5A) + 5{z - Xa){Pa,ga''^9a) + Kz - Xbj^b^Qb^'^gb) + S{z - Xc){Pc, 9c^ Sgc)], 

c+a,h 

(29) 

where <> is the Killing form, b is the exterior differentiation operator, 
-A and A are holomorphic and antiholomorphic parts of the connection cor- 
respondingly, 

-{pa,ga) £ T*Ga where T*Ga is the cotangent bundles at the marked point, 
Pa G Lie*{Ga),ga ^ Ga- Symplectic form is invariant with respect to the 
action of the group Gi: 

G, = {f{z,z)eC^{j:),G} (30) 
The corresponding gauge transformations have the following form: 

A = fLf-' + fdf-\ 

A = flr^ + fdf-' (31) 

Pa — faPafa > 9a — 9afa ■ 



The moment map equation takes the following form: 



dL = 5{z- Xa)Pa + S{Z- Xb)pb + ^ 6{z - Xc)p, 



cf=a,o 



After transformation (16), we get: 



dL = d{z - Xa)p° + 5'{Z- Xa)p^ + ^ S{Z - Xc)Pc- (-33-) 

c^a,b 

Solving this equation, we get the Lax operator with the second order pole. 
Then let us find such transformation of the form which gives us the new 
form ^new corresponding to (33). Let us introduce the logarithmic coordinates 
In g on the group Sl{N, C), which are the coordinates on the algebra sl{N, C), 
then In^a = Xa, lngh = Xj, and Qa^Sga = SXa, g^^Sg^ = SX^, and consider 
the following transformations: 

Xa = X°e-^ - X\ Xb = -X^e-^ - X\ (34) 

Let us define the coupling between the algebra X and the coalgebra P as 
Res-iTr{PX). Substituting (16), putting ao = 1) ai = —1, /3o = 0, A = 
1, and (34) to the form after expansion in e, collecting similar terms and 
returning to X^ = {g^)~^5g^, we get the new form ^new'- 

n= [ [{dAAdA)+5{z-Xa){p^,go^dgQ) + S'{z-Xa)(j)\go^5go)]+ V (pc,5c" 

(35) 

It is possible to get (33) from this form. 



3 Systems arising from the elliptic Gaudin 
system 

Hamiltonians of Gaudin system in the case of two point fusion 

Let us begin (as in the rational case) from the two point fusion and consider 

two marked points Xa and X5. According to what was said in the introduction, 
the integrable system is defined on the symplcctic factor-space R ([4]) which 
has the following form in the case of two points: 

R={u,v)x{OaxOb//H), (36) 

where Oa and are the coadjoint orbits of the group SL{N, C), 

H is the Cartan subgroup, 

u,v ^TL belong to Cartan subalgebra. 

The moment constraint corresponding to the action of the remnant gauge 
group has the following view: 

Ep'' = (37) 

a 

Let us consider the case N = 2 (i.e. p G sl{2,C)) and get the Hamiltonians 
by decomposing < > into the sum of Eisenstein functions, where L is 
the matrix of the Lax operator in the holomorphic representation, with the 
diagonal part 

L'' = v' + J2p^-E^iz-Xa) (38) 

a 



and non-diagonal part 

^'^' = EPa •'^(«'^^-a^a), (39) 

a 

where: 

/ i-i ^ e{u'^ + Z - Xa)e' (0) 

0(u^^)0[Z — Xa) 

In the case of the two point fusion we get: 

< >= (v^ +pI^ ■ Ei{z - Xa) +pI^ ■ Ei{z - Xb)f+ 
+pf . Ei{z - Xa) +pf ■ Ei{z - Xb))^+ 
+plY^\-E[iz - Xa) + e[{u'^)) +pfpl\-E[{z - X,) + E[{u'^))+ 

+'^V^aPb' • ^{U^'^^ Z - Xa)^{u^^,Z - Xb) + 2pl^pl^ ■ ip{u^^ , Z - Xa)^{v}'^ , Z - Xb). 

(41) 

Analyzing zeroes and poles in < > we get the decomposition in the fol- 
lowing form (Appendix C): 

< >= E(^i(^ - • + El{z- Xa) ■ H^a'') + H^a'° (42) 

a 

, where: 

h!''=<pI>, (43) 

H^'' = v'Pa + vVa' + \PaPl' + pfpf] ' E^{Xa " Xb) + 

|0 12„21 ^(^^^+^fe-^a)^'(0) , .J1112 0{u^'+Xa-Xb)e'{Q) (^4) 
^PaPb ^(^12)^(^^_^^) ^ PaPb ^(^12)^(^^_^J ' 

H^'^ = + - IpI'pV +Pfpf] ■ {E[{Xa - Xb) + Elixa - Xb))+ 

+ lplVa'+plVb']-E[{u'')) + 

^2p-piw) - E.iu- + - .„)) . 

^2pfpiw) - E^iu- + xa - Xb)) . 

(45) 

Let us calculate the phase space dimension of the integrable system. There 
arc eight variables: u^'^ , , p^^ , p^ , p"^^ , pl^ , pl"^ , pf^ in the system which is in 
accordance with two Casimir operators and three Hamiltonians. Fixing the 
value of Casimir operators we choose the orbits i?^'^ = (Aa)^, H^'^ = (A;,)^ 
and taking into account (37) we get: 

/ 11\2 , 12 21 /■> \2 / 11\2 , 12 21 /-> \2 11 , 11 n /'/lft^ 

[Pa ) +PaPa = (^a) , [Pb ) + Pb Pb = (Afc) , Pa + Pb = ^ (46) 

Note, that the remnant gauge group in the general case (not necessarily in 
the holomorphic representation) consists of the doubly periodic functions of 
(z, z) on the diagonal because they do not change the gauge fixing in the mo- 
ment map equation (sec Symmetries). Taking the Fourier-series expansion of 
these functions we get the basis in the space of the remnant gauge transfor- 
mation. After that let us consider the gauge fixing being in accordance with 



the coadjoint action of the remnant gauge group on p whieh are independent 
of (z, z) diagonal matrices being the Cartan subgroup with the corresponding 
moment constraint (37). So we must consider functions which have only zero 
harmonic cq = exp(+Q;) in the expansion and generate one-parameter class. 
The invariant with respect to the gauge fixing variables take the form: 

12 21 _ 12 21 _ 12 21 _ 12 21 _ 11 11 

Pa Pa Ph Pb =y, PaPb Pb Pa =^2, Pa , Pb , U, V 

(47) 

Finally, we get: 

{pl'f + x = {Xaf, {pVf + y = {Xb?, pI'+pV=0, xy = ziZ2 

(48) 

It is clear that six variables (all except u, v) can be expressed through two ones, 
for example zi,Z2- So, we get four-dimensional phase space in accordance with 
two Hamiltonians. Note, that the Hamiltonians, as functions on this phase 
space, are invariant with respect to the action of the remnant gauge subgroup 
because phase space is the factor-space with respect to the action of this 
subgroup. 

Hamiltonians of the system obtained by the two point fusion 

Let us realize the point fusion in accordance with the transformation (16). 
Note, that commutation relations between (jpy^ and (p^Y^ do not depend on 
the elliptic function and are the same as in the rational case. The additional 
issue in the evaluation of these limits is the decomposition of the elliptic 
function into series of e. In the case of p G sl{2,C) after calculation, setting 
= u and = —v^ = v, we get (Appendix D.): 

^ < Ll^^ >= E^{Z - Xa) ■ H^'^ + E';{z - Xa) ■ i?"'" + E[{z - Xa) ■ H^'^ + H''^\ 

(49) 

Where : 

= -a + (c - 6) • Ei(ti) + ^ • [e[{u) + Ei{uf] - 2vg, 
//2,o = y2^g^. [^e';'{u) + iE[{u)f] + a-E[{u) + ib-c)- [^E'liu) + E^iu) ■ E,{u)]- 
• [^<'(«) + iE[iu)f + e';{u) ■ Eiiu) + E[iu) ■ Eiiu)% 

Qj O 

(50) 

where we use the following notation for the variables invariant with respect 
to the gauge fixing: 

{pr\pr'=a, {py\p'?'=h, {p'r{p'r=c, {p'r{pT=d=-, {p'r=9 

(51) 

Let us calculate the dimension of the phase space of the new system. There 
are eight variables in the space {u,v) x {p^Y^ ■ .ff^'^ and H^'^ are the Casimir 
operators. The moment map for the remnant gauge action takes the form 
(pO)ii = 0. So, we get the four-dimensional phase space and two Hamiltonians. 
Fixing the Casimir operators gives the equations: 



\- g = Ao, b + c = Xi 

a 



(52) 



So we get four dimensional phase space and two Hamiltonians. The Poisson 
brackets in the new variables take the following form: 



{v,u} = l, {c,b} = 0, {a,g} = c-b, 

he he 

{c,a} = -2ag, {b,a} = 2ag, {b,g} = , {c,g} = —. 

a a 

The Poisson brackets for cZ = ^ take the following form: 



(53) 



be 



rbc , , rbc , rbc bc rbc . 

{-,a} = 2g{b-c), {-,3} = 0, {-,b} = -2g-, {-,c} = 2g- 
a a a a a a 

(54) 

Symmetries Hamiltonians defined on the phase space must be invariant 
with respect to the remnant gauge action - the Bernstein-Schvartsman group, 
which is the semi-direct product of the Weyl group and the lattice shifts. 

Let us consider (in accordance with [4] ) the gauge transformations preserv- 
ing the gauge fixing L = diag{sl[2, C]), chosen in the moment map equation: 



dL + 



27ri 



• [L, L] = 2m ■ {6^{xa)pa + 5'^{xb)Ph 



T — T 

with boundary conditions 

L{z + l) = L{z), L{z + t) = L{z), 

where: 

L = 



u 




u 



(55) 

(56) 

(57) 



Then the remnant gauge group in the general ease consists of the doubly 
periodic Cartan valued functions of (z, z) because they preserve the gauge 
fixing in the moment map equation. Taking the Fourier-series expansion of 
these functions we get basis in the space of the remnant gauge transformation 
which consist of harmonic having the following form: 



f = exp 



„ . / iZ - Z jTZ - TZ, 

2m ■ (m* + n* —] 



T — T T — T 

Let us consider the transformation of L under the basis function action: 



(58) 



2m- 



1 



r — r 



u 



f ■ 2m- 



1 



r — r 



u 



n + fdf 



-1 



where: 



/ 



f 
f 

Prom this we get law of the transformation for n': 

^ + m* - n V 
Let us consider the transformation of u^"^ in the following form n^^ 



(59) 



(60) 



(61) 



u + 1. We may represent this transformation as the following two: 



u 



^ ^ and — > v!^. Then / take the form {m^ = l,n} = 0): 



/ 







, = exp 



2m 



z — z 



T — T 



(62) 



It is possible to define the action of / on p^. 



o ^ faPafa 



fa J a 

We obtain the transformation law for p^J^ : 



22 
Pa 



(63) 



Pa^Pa^ Pl^^Pa^W 



T — T 



Pf Pf exp 



-2m 



T ~ T 

(64) 

It is possible to obtain similar expressions for v)- ^ v)- + r, in the general case 
we have: 



P^'i p'i ■ exp 



u^^ + - nV 



p^: - pz 



(65) 



r — T r — r 

ip{v}^ ,z — Xa) — > ^{u^'' , z — Xa) ■ exp +27rm* • {z — Xg) 



Transformations of the system obtained by point fusion Let us 

do the gauge transformations oi L,L,pl^ in the moment map equation using 
the following function: 



/ 



exp 




—2Triu 



z — z 







exp 



-27riu 



z — z 



T — T 



(66) 



So we get the moment map equation in the following form: 

BL = 27ri • {5'^{Xa)Pa + 5^{xh)ph), 
with the boundary conditions for the matrix of the Lax operator: 



{z+l) = {z) , {z + t) = (z) ■ exp 



-2mu^^ 



(67) 



(68) 



Nondiagonal part of the matrix of the Lax operator take the holomorphic 
form: 



^" e{uij)e{z - Xa) 



e{u^^ + z - Xb)e' {Q) 

e{u^^)9{z - Xb) 



(69) 



It is necessary to express the new variables through old ones {pa)new — 
Pa ' r\f^^)~^ to find the transformations for the new variables {paTJew 
accordance to the transformations n*-^ n*-^ + — nV. Having found the 
eventual transformation we get it in the holomorphic gauge: 



u^^ + m* - nV 



pV 



P's ■ exp 



-27rz ■ n'x, 



(70) 



Now we can consider the transformations of the system obtained by point 
fusion. Having done the over scaling transformations of p^^ in the matrix of 
the Lax operator in the holomorphic gauge and at the same time the trans- 
formations of the variables under the shift u^^ , we get the modified action of 



the Bernstein-Schvartsman group on the new variables {p^Y-' , (p{u^^,z — x^) 
and its derivative, having the following form: 

u^^ = u ^ u + m — riT, 

(p{u, z-Xa) ^ (p{u, z - Xa) ■ exp [+27rm • {z - Xa)] , 

if {u,z — Xa) — > [— 27rm ■ ip{u, z — Xg) + f {u, z — Xa)] ■ exp [+27rm • {z — Xa)] 

(71) 

Matrix of the Laix operator of the new system Let us consider 
the matrix of the Lax operator in the holomorphic gauge in the case of two 
points . Having done transformations (16) and the decomposition in e of the 
resulting expression, passing to the limit e — > we get: 
diagonal part - 

L'' = v'-{p'f-E[{z-Xa) (72) 

nondiagonal part - 

Vi = (jPyi . <^(u*^-, z-xa)- {p'f ■ ^k^'\ ^ - ^a) (73) 

The matrix of the Lax operator and the corresponding Hamiltonians depend 
on the method of the decomposition of p*". It means that it is possible to con- 
sider decomposition in e and e. In the present work we consider holomorphic 
decomposition parameterized by e. 

In the elliptic case with two marked points the total moduli space of the 
system is a fibered space with the base defined by the elliptic module r. The 
fiber of this space is a torus with marked point. This marked point corresponds 
to a new system having there a pole of second order in the matrix of Lax 
operator. This follows from the considering the module xy, — Xa (with fixed 
r). In the other words we can fix, for example Xa and then X5 will take value 
in the torus. When Xb — > Xa the new system obtained by the fusion of Xh and 
Xa corresponds to the fixed point. 

Generalizing to the case of r point fusion we get the matrix of the Lax 
operator of r point fusion in the following form: 

diagonal part - 

= v' + (pi)" . k2 ■ e[{z - Xa) + ... + {pT ■ K ■ E^{\z - Xa) (74) 
nondiagonal part - 

L'^ = (jPyi . ~ki ■ ip{u'^,Z -Xa) + ... + {pn'' ■ ~kr ■ ^^'^ {u'^, Z - Xa), (75) 

where kr is similar to fc^, commutation relations are the same as in the 
rational case. Note, that the symplectic structure is the same as in the rational 
case, too. Integration in the formula (35) is taken over the torus. 

4 Inozemtsev Limit 

In the Inozemtsev Limit procedure we must satisfy the requirement of the 
absence of singularities in the integrals of motion which are the traces of the 
Lax operators which coincide with the solution of the moment map equation: 

\<L^>,\<L^>, ■■■,1<L'> (76) 



Taking the Eisenstein-series expansion we get required number of Hamiltoni- 
ans as the coefficients of the decomposition. Let us write the first nonvanish- 
ing term of the Eisenstein functions in the Umit UI2 ^ oo, taking into account 
u = u + tL02, where t is some parameter: 



function 


t = 


< t < 1 


t = 1 


l<t<2 


Ei{u) 


icosh(f) 
"^2 smh(f) 


+ \ 


+ \ 


+\ 


E[{u) 


1 1 
4smh2(2) 


_^-U-tui2 


-2cosh(ii) • 


_^+u-{2-t)uJ2 


E'liu) 


1 cosh(f ) 


_j_g-M-tl^2 


-2sinh(n) • e"'^^ 


_^+u-{2-t)u)2 


E';'{u) 


11 3 1 
4sinh2(f) 8sinh4(|) 


_^-U-tul2 


-2cosh(u) • e-'^2 


_g+w-(2-t)a)2 












e[''\u) 


1 cosh(f ) 
''"4sinh3(|) ■•■ 


_j_g-U-tu)2 


-2sinh('u) • e^'^^ 


_g+u-(2-t)a)2 


E^''+'\u) 


1 1 

4sinh2(f) - 


_^-U-tu>2 


-2cosh({t) • e~'^2 


_g+W-(2-t)w2 



Here we put: 



-in, Im{u2) = 

Degeneration of Ei{u) up to the second unvanishing order has the form: 



(77) 
(78) 



1 



= Yl coth(- -A:w: 



0<t<l: ^ + ^e-^-^^ 
t = 1 : ^ — sinh('u) • e"'^^ 



(79) 



l<t<2: 



1 1 



,+u-(2-t)a;2 



According to this data we will consider the degeneration of and H^'^. 
First of all let us consider some general remarks. We take into account only the 
first term in the decomposition Ei (u) because the next terms must vanish in 
order to avoid the appearance of singularities. So if < t < 2 we get Ei{u) = ^ 
after the ovcrscaling. The terms in H^'"^ containing E'i{u), containing (£'^(u))^ 
in if^'^, and 2vg in H^'^ also vanish up to requirement of the absence of 
singularities after overscaling H^'^. The variables a,b,c,g and z transform in 
the following way: 



e-x.u;2^^ 



(80) 



^^^-Xau>2^^ ^^g-X6a;2^^ g^g-Xca;2g^ g 

Z Z + SUJ2 

There are the following possible relations between t and some parameter s: 

s = t,2-t (81) 

Let us consider the case of preserving all Lie-Poisson brackets between a, b, c, g 
after overscaling. The conditions of the absence of singularities in the brackets 
after overscaling have the following form: 



Xb > Xg, Xc > Xg, Xg > Xb - Xa, Xg>Xc- Xa 



(82) 



The equality sign corresponds to existance of the Lie-Poisson bracket limit, 
so: 

Xb = Xc = Xg, Xa = (83) 



And the brackets take the following form: 

{v,u} = l, {c,b} = 0, {a,g} = c-b, 

be be (84) 

{c, a} = -2ag, {b, a} = +2ag, {b,g} = - — , {c, 9} = — 



Casimir operators take the following form: 



ij2'^ = — + 5^ H^'^ = b + c, 



be 



2,3 



(85) 



and Hamiltonians in accordance to the value t take the following form: 



H"^,^ = {c-b)-^ + — -^, 0<t<2 



1 be 1 



H 



2,0 



.-.I be 



2 a 2' 
1 



V — g ■ cosh({t 



6 a 
1 



12' 

c).e- + -.[-e^ 



12 



e-% t = 1 



v'^-f. e+«^ - (6 - c) . e+" + - ■ e+"^, 1< t < 2 



(86) 



5 Appendix A. 

There are main formulas for elliptic functions in this section (they are bor- 
rowed from [4] and [9]). First of all we define ??-function: 



i?(z,r) = gie-f(e^'^^ - e"''^^) - g'')(l - g"e2'^^^)(l - ^"e-^'^'^), q = e^""^^ 

n=l 

(87) 

where r is the complex module, and Eisenstien functions: 
Ei{z,t) = dJog'&{z,T), Ei{z,t) ^ + ... 

E2{z,t) = -d,Ei{z,T), E2{z,t) « 4 + ... 



E!,iu) 



Ei{u + v) + Ei{u-v)- 2Ei{u) 



E2{u)-E2{v) 

Their relations to the Weirstrass functions have the following form: 

az,T)=Ei{z,T)+2rii{T)z, p{z,T)=E2{z,T)-2rii{T). (89) 

where 771 (r) = C{^)i ^^id ({z) is Riman ^-function. Eisenstien functions have 
the following representations: 



Ei{z) = i^ E coth(J-fca;2), E2iz) = - ^ 



k=—oo 



=-jy: 



^fc=_^sinh^(i-A;a;2) 
cosh ( I — kuj2) 



^^^^smh\^-ku;2) 



(90) 
(91) 



There is the following expression in the matrix of the Lax operator: 



It has a pole at the point 2 = and res|^=o</^(u, z) = 1. 
The important identity connected with ip: 

(p{u,v)ip{-u,v) = E2{v) - E2{u), (p'^{u,v) = (p{u,v){Ei{u + v) - Ei{u)) 

(93 

Behavior on the lattice: 



(94) 



e{z + 1) = -e{z), e{z + r) = -q-^2e-^^'^0{z), 
Ex{z + 1) = Ei{z), Ei{z + t)= Ei{z) - 27ri, 
E2{z + l) = E2{z), E2{z + t) = E2{z), 
(p{u + l,z) = (p{u, z), (p{u + T,z) = e~'^'^^^(p{u, z). 

6 Appendix B. 

Let us consider the canonical symplectic form u on T*G = G* x G { G E 
SLiN,C)y. 

LJ = 5<p,g ^5g> (95) 

where <> denotes the Kilhng form on Lie(G), p € sL*{N, C) g £ SL{N, C), d 
- external differential operator. We will find the commutation relation between 
F{p,g) and H{p,g) defined on T*G. There is Lie-Poisson bracket on T*G 
([11]): 

{F,H} = ctx'djFdkH, (96) 

where C/*^ is the structure constants of the Lie algebra Q 
dj = g|j x' is the coordinate in the space Q*. 

We want to rewrite (96) in more convenient form for the next calculations 
taking into account the explicit dependence of F and H oi p and g. Let us 
represent (96) in the following form: 

{F,H} = XfH, (97) 

where Xp is the vector field corresponding to F, so: 

arp air 

{F, H}=<P,[^,-Q^]>-<9, {F, H} >, (98) 

where {F,H} = ^^-^^. Put H=<p- Eji > and F =< p ■ Eik >, so 
we get the formula (1): 

{F,H}=pkjS''-PiiS''^ (99) 
In the case of p G sl{2, C) and two point fusion: 



G = 



p\ 



(100) 



Using (98), we get 

r 0,27 Ijfc/i cil k,i+N xki i,l+N 

iPi ,Pi } = S Pi'' -S 'p{ 
1 (101) 



7 Appendix C. 



Decomposition of < L >. Put L in the following form: 



L 

then: 



^21 ^22 



(102) 



iTrL2 = (Lii)2 + Li2^2i (103) 

Let us consider the decomposition of the diagonal and nondiagonal parts of 
< > into the Eisenstien functions series in the case of two point. 

Diagonal pcirt In this part we will use the property of elliptic functions 
with equal periods to differ in constant term if they have equal poles with the 
main parts in the parallelogram of periods. 

Lii = + ^P"^ii ■Ei{z- Xa) + {pb)ii ■Ei{z- Xb)] (104) 

i 

Let us write the supposed decomposition for TrLf^ (we do not write i, sup- 
posing summation): 

v'^+pI- Ef{z - Xa) +pI ■ Ef{z - Xb) + 2vpa ■ Ei{z - Xa) + 2vpb ■ Ei{z - Xb) + 

+2paPb ■Ei{z- Xa)Ei{z - Xb) = 

= -pI ■E[{z- Xa) -pl-E[{z-Xb) + a- Ei{z - Xa) + (3 ■ Ei{z - Xb) + cq + v^, 

(105) 

where: 

a = 2paPb ■ El {Xa - Xb) + 2vpa, (3 = 2paPb ■ Ei {xb - Xa) + 2vpb, (106) 

Now we must define cq. In order to do this let us consider z ^ Xa, representing 

z = Xa + e, so as e ^ 0. Expressing cq from (84), decomposing all terms in e 
series, collecting the similar terms and using the condition {pa)ii + (pb)^ = 0, 
we get: 

Co = +PaPb ■ [E'liXa - Xb) + El{Xa - Xb)] (107) 

Note that in our calculation cq does not depend on the chosen point (on the 
way of calculation). 

Nondiagonal part In this part we will use the property of elliptic function 
with equal periods to differ only by constant multiplier if they have equal zeros 
and equal multiplicity of the poles in the parallelogram of periods. 

The term of nondiagonal part of < > in the case of two point is 
determined in formula (41) of chapter 3. These are the following expressions 
having f?-functions: 

(p{u,Z - Xa)^{-U,Z - Xa), (p{u, Z - Xa)<fi-U, Z - Xb) , (108) 

where: 

f ij ^ 0{u'^ + z - xs)e' jo) 

(flu-' ,Z — Xg) = r (J-Uyj 

and s = a,b. Analyzing zeroes and poles of this expressions we will extract the 
elliptic parts, so this enables us to obtain the decomposition in the Eisenstien 
functions series depending on z — Xg- 



Decomposition (p{u^ z — Xa)<f{—u, z — Xg) 

9{U + Z- Xa)9{-U + Z- Xa){e'{0))'' , . 
(p{u, Z - Xa)(p{-U, Z - Xa) = , , , rr^ (HOj 

Zeros: 

U^±{z-Xa) (111) 

Poles: 



■U ^ 0, i{}{u,Z - Xa)i~p{-U,Z - Xa) ^ \ 

Z^Xa, (p{u, Z - Xa)(p{-U, Z - Xa) ^ + 



(112) 



{Z - Xa)'^ 

It follows from the decomposition of ^-functions with respect to small param- 
eter. Thi decomposition takes the following form: 

ip{u, Z - Xa)ip{-U, Z-Xa) = E[{u) - E[{z - Xa) (113) 

The similar expression we get for z — Xf,: 

(p{u, z - Xb)ip{-u, z-Xb) = E[{u) - e[{z - Xb) (114) 

Decomposition <^(u, z — Xa)ip{—u, z — Xb) 

9{u + z-Xa)9{-u + z-Xb){9'iO))'^ . . 

^{u,Z-Xa)^{-U,Z-Xb) = (115) 

Zeros: 

U^-{Z-Xa), U^+(Z-Xb) (116) 

Poles: 



n — > 0, f{u, z — Xa)(p{—u, Z — Xb) ~ 



u 



Z^Xa, (p{u,Z - Xa)(p{-U,Z - Xb) +C ■ (117) 

[Z -Xa) ^ ^ 

Z ^ Xb, (p{u,Z - Xa)(p{-U,Z - Xb) ^ -C 



{Z - Xb) 

Using these zeros and poles we may define the combination of Eisenstien 
functions: 

(p{u, Z - Xa)(p{-U, Z - Xb) = C ■ [-El {z - Xa) + Ei{z - Xb) - Ei (u) + Ei{u - Xa + Xb)] 

(118) 

We may define the coefficient C comparing coefficients in the left and right 
parts, for example at z — Xb'- 

(J ^ 9iu + Xb-Xa)9'{0) ^^^g^ 
9{u)9{Xa - Xb) 

Finally we get: 

(p{u, Z - Xa)(pi-U, Z - Xb) = 

= ^^'1,^'!'"'""^^^°^ • [-Ei{z - Xa) + Ei{z - Xb) - Ei{u) + Ei(n - x„ + Xb)] 
9[u)9{xa - Xb) 

(120) 

and for (/p(n, z — Xb)<p{—u, z — Xa)- 

(p{u, Z - Xb)ip{-U, Z-Xa) = 

= • - ^^b) + El{z - Xa) - Ei{u) +Ei{u-Xb + Xa)] 

9{u)9{Xb - Xa) 

(121) 



8 Appendix D. 

Obtaining Hamiltonians of the new system. Let us find: 



^<L'>={L^^f + L^''L'\ (122) 



(124) 



in the case of the two point fusion. Let us consider L^^L?'^: 

= [(/)12 . ^(^^ z-Xa)- ip')'' ■ <p'{u, Z - Xa)] ■ [(p°)'^ • <p{-U, Z - Xa) - {p')'' ■ f' {-U, Z - Xa)] 

(123) 

We need the next formulas to calculate (123): 

ip{+U, Z - Xa)^p{-U, Z - Xa) = E[{u) - E[{z - Xa) 
ip' {+U, Z - Xa)fp' {-U, Z - Xa) = 

= -[{e[{u))'' + EI{u)e[{u) + E,{u)E';iu) + 

+ [Ef{u) + E[{u)] ■ E[{z - Xa) - ^E';'{z - Xa) 
ipi+U, Z - Xa)'p'{-U, Z -Xa) = 

= -[Ei{u)e[{u) + ^E'liu)] + Ei{u) -E'liz- x„) - ^e';{z - Xa) 

(p{-U, Z - Xa)(p' i+U, Z - Xa) = 

= +[Ei{u)e[{u) + ^E^iu)] - Ei{u) -E'liz- x„) + ^e'^^z - x„) 

Put them in (123). Having done the rearrangement and summing with 
the decomposition of (Ln)^ + (-^^22)^, we get: 

^TrL' = E'^iz - Xa) ■ H''' + eI{z - x„) • H^'^ + e[{z - Xa) ■ H^'' + H''' 

(125) 

Where : 

H''' = -1{p')'Hp')'' - \{p'f\p')'\ H^'' = ^iiprW + iP^'ip')''] 

(126) 

^2,2 ^ _(p0)12^^0)21 + . [_(^0)12(^1)21 + (^0^21(^1^12] + 

+ip')'\p')'' ■ [E[{u) + Efiu)] - 2v{p'f^ ^^^^^ 



= v^ + UpT)' ■ [Ie: in) + iE[iu))^] + {py\py' ■ E[{u)+ 
+[(p'f\p')'' - ip')''ip')''] ■ [Ei{n)E[{u) + ^<(«)]- (128) 
-{p')'Hp')'' ■ [{E[{u))' + EI{u)e[{u) + E,{u)E';{u) + ^E'^iu)] 
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